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Abstract 

We present a straightforward model of cosmic-ray propagation in the Galaxy that can account for 
the observed cosmic-ray positrons entirely as secondary products of cosmic-ray nucleons interacting 
with the interstellar medium. In addition to accounting for the observed energy dependence of 
the ratio of positrons to total electrons, this model can accommodate both the observed energy 
dependence of secondary to primary nuclei, like Boron/Carbon, and the observed bounds on the 
anisotropy of cosmic rays. This model also predicts the energy dependence of the positron fraction 
at energies higher than have been measured to date, with the ratio rising to ~0.7 at very high 
energies. We briefly point out the differences between this model and the model currently in 
wide use that does not account for the observed positrons as secondaries and so prompts the 
interpretation of the observations as evidence for alternate origins of positrons. 



INTRODUCTION AND OVERVIEW 

The recent measurement of the positron fraction R e+ (E) = F e+ /(F e+ + F e J) at energies 
E up to 300 GeV by the AMS collaboration p] is an important contribution to cosmic-ray 
physics and poses a challenge to predict F e+ (E) with similar precision. 

A striking feature of the AMS data confirming, with high statistics, earlier observations 
[21 [3] with unprecedented accuracy, is the monotonic increase of R e+ (E) from ~0.052 at ~10 
GeV to ~0.155 at ~300 GeV. This contradicts an earlier calculation of R e+ (E) expected 
assuming positrons are secondaries produced by primary cosmic-ray nuclei in the interstellar 
medium, by Moskalenko and Strong (MS) [4J which predicted a monotonic decrease (see the 
lower panel in Fig. 1). Analysis with the MS model seems to require new physics to explain 
the discrepancy. It has been proposed that the positron excess could either originate from 
pulsars' magnetospheres [5] or from the annihilation or decay of dark matter |Bj, a more 
exciting explanation which is constrained by the absence of high energy gamma rays from the 
center of the Galaxy (?]. A kinematical cut-off in the positron spectrum at E ~ (| — |)M(x), 
with M(x) being the mass of the dark matter, is generally viewed as an indication for dark 
matter. 

We show here that the observed decrease of R e+ (E) up to ~6 GeV and its subsequent 
increase can be explained entirely as cosmic-ray secondaries, if positrons have a 1 — 2 Myr 
residence time in the general interstellar medium, independent of their energies. Our cal- 
culations predict a dip in the positron fraction beyond 300 GeV followed by an increase 
beyond 1000 GeV to reach an asymptotic value of ~0.7. This prediction is shown in the 
lower panel of Fig. 1. The difference of up to a factor ~10 at 300 GeV between our cal- 
culated F e+ (E) and that predicted by MS [4 J indicates that present theoretical calculations 
are model-dependent and do not allow us to reach firm conclusions [8HT3]. including that 
decay or annihilation of massive dark matter in our Galaxy is responsible for the positron 
spectrum. 

The rest of this letter is devoted to a presentation of the analysis leading to our results and 
an explanation of the differences between the models for cosmic-ray propagation suggested by 
the positron observations and those currently in vogue for analyzing cosmic-ray observations. 
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FIG. 1. Upper panel: The solid black line represents our fit, Ft(E), to the spectrum of the total 
electronic component observed in cosmic rays, the dotted line shows the source function S e +(E), 
with riH = 1 cm" 3 and r = 1 Myr, which fortuitously lies very close to the data points representing 
the observed positron spectrum Fams obtained by multiplying AMS-02 data on positron fraction 
by Ft(E). The dashed line represents F e+ (E) the theoretical spectrum with S e +(E) as the source 
function, including propagation effects during a residence time r = 2 Myr and n# = 0.5 cm" 3 . 
Lower panel: Our predicted positron fraction, R e+ (E) = F e+ /Fx, with uncertainties is shown; the 
shaded steeply falling region is due to MS models [J| . 

ANALYSIS OF THE MODELS. 

The spectral intensities of cosmic-ray positrons provide a hitherto unavailable probe for 
the study of the origin and propagation of cosmic rays. First, they are not ubiquitous like 
electrons, which may be accelerated to cosmic-ray energies in the sources; positrons have 
to be generated either as secondaries through the n + — > ji + — > e + + v e + and other 
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decay chains of mesons produced in the cosmic-ray interactions in the interstellar medium 
or should be generated in exotic processes in pulsar magnetospheres or through the decay 
or annihilation of dark matter. Secondly, in high energy cosmic-ray proton interactions, the 
positrons carry away a small fraction ~3-6% of the energy of the primary, unlike secondary 
nuclei, like B, which emerge with the same energy per nuclei as their parent nuclei, C, O, 
etc. Thus the observed spectrum of positrons will carry a signature of its origins. 

We have calculated the spectrum of positrons, F e+ (E), by multiplying the positron frac- 
tion, R e+ (E), measured by the AMS instrument by a fit to the spectrum of the total elec- 
tronic component, Ft(E). Both these spectra are displayed in Fig. 1. The PAMELA results 
[2] are consistent with F e+ (E) displayed here. The spectrum of positrons F e+ (E) has the 
form A ■ E~P+ with /3 + ~ 2.65, almost identical with that of the total nuclear component of 
primary cosmic rays. In contrast, the total electronic component has a spectrum that has a 
spectral index (3 T « 2.2 below a few GeV, steepening to an index of (3j> ~ 3.1 until ~1000 
GeV, beyond which there is a rapid decrease of the intensities. 

The rate of generation of positrons in the interstellar medium depends on the well- 
measured intensities of the nuclear component of primary cosmic rays and the cross sections 
for meson production, which exhibit a scaling behavior with adequate accuracy. A recent 
compilation of cosmic-ray proton and He spectra at high energies by Benard et al. [14), 
yields a hard spectrum for the nucleons with spectral index ~2.65 in the relevant energy 
region. Based on this spectrum and the calculations of MS (see Fig. 4 in |3j), we estimate 
the source spectrum of the positrons to be 



with So ~ 4x 10 _4 cm~ 2 s~ 1 sr~ 1 GeV _1 Myr~ 1 and (3 + ~ 2.65, when n# is given in terms of the 
number of hydrogen atoms cm -3 and E is in GeV. We display in Fig. 1, S e+ (E) for n# = 1 
cm -3 , which is a remarkably close fit to the observed spectrum in the energy band 3 — 100 
GeV, with a fortuitously close normalization. With n# = 1 cm -3 , 1 Myr corresponds to 
an effective grammage, A e+ w 1.7 g-cm -2 for the positrons. In order to obtain the steady 
state spectrum of positrons, we should account for the energy loss suffered by positrons 
through synchrotron radiation and inverse Compton scattering against the 2.7 K microwave 
background during their residence time in the Galaxy. This loss is parametrized as 



S e+ (E) = S n H E-P+ 



(1) 



dE 
~dt 



-bE 2 



(2) 
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where b pa 1.6x10 3 GeV 1 Myr 1 and E is in GeV [9j. The steady state spectrum of 
positrons [TS] in the interstellar medium, F e+IS M(E), is 

• 1/bE 



e+ISM 



(E) = / S e+ (E) ■ (1 - bEt)^+- 2 e t/T dt 
Jo 

pa 4 x 10- A n H rE~ 2m for 1 GeV < E < 300 GeV 
4 x 10~ 4 



1.656 



n H E~ 6 - b5 for E > 300 GeV. (3) 



Here r is the leakage lifetime for the positrons in Myr units. As positrons enter the he- 
liosphere and propagate to near-Earth locations, they suffer solar modulation. Accounting 
approximately for these effects, which affect only the very lowest end of their spectrum, we 
get 

F e+ (E) = F e+ISM (E)e- E -/ E (4) 

and display this spectrum in Fig. 1 for n#= 0.5 cm -3 , r = 2 Myr and E m = 0.5 GeV, 
which provide a good fit to the data. The product n#r pa 1 corresponds to a grammage 
A e+ pa cm#n#r pa 1.7 g-cm -2 . Note that the energy losses have steepened the high energy 
part of the positron spectrum to E~ 3,65 and has made the spectral intensity independent of 
the leakage lifetime, a result that can be shown to be true even if r were to be dependent on 
energy. Note that the calculated spectrum of secondary positrons fits the observations well. 
The free parameters in this fit are the leakage lifetime, r pa 2 Myr and the mean interstellar 
density n H = 0.5 cm -3 . In the lower panel, the positron fraction, R e+ (E) is displayed and 
fits the observations, as expected, because the e + spectrum agrees well with the calculation. 

What is the expected behavior of the positron fraction at higher energies? The answer 
to this question rests on the following considerations: In the energy region up to ~2000 
GeV, we have the observations of the total electronic component. Accordingly, the ratio 
R e+ (E) = F e+ (E) / Ft(E) is easily calculated. The sharp steepening of Ft(E) beyond 1000 
GeV is attributed to the discrete nature of the cosmic-ray sources and the energy losses 
suffered by the electronic component in the finite amount of time needed for them to arrive 
at the Earth [SJ El CE2J US! ED]- Even though the primary electronic component cuts off, 
the secondary electrons and positrons that are produced in the interstellar medium that 
surrounds the Earth suffer only the aforementioned steepening and continue as ~^~ 3 - 65 
at least up to 10 4 GeV. At these energies they dominate the flux and their ratio will be 
controlled by their production characteristics. Keeping in mind that there exists an excess 



of protons over neutrons (bound in nuclei) in the primary cosmic-ray beam and the fact 
that inelastic diffraction in the forward direction dominates the secondary cosmic-ray flux, 
we may note that e + is favored over e~ in the production and the secondary e~ is about 
~0.5 of the e + . Thus we expect the R e +(E) to reach ~0.7 at E > 10 3 GeV. 
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FIG. 2. The grammage A(E) for the various cosmic-ray particles as a function of the kinetic energy 
per nucleon or per positron is displayed for the leaky box model of Davis et al. [T7| as dashed lines 
and estimated from the diffusion model of MS [4J as solid lines. In our model, much of the energy 
dependent part of the grammage is attributed to traversal in the sources and a constant value ~ 
1.7 g-cm~ 2 , independent of energy, is traversed in the interstellar medium of the Galaxy. 

DISCUSSION. 

It would be appropriate now to address the question: what are the significant differences 
between the two models, one by MS [I] and the other described here, both displayed in the 
lower panel of Fig. 1, that they make such diverse predictions for the positron fraction? 
The motivations for the cosmic-ray modeling has been provided by the observations of the 
ratio of secondary nuclei like B to that of their parent nuclei like C and O. Once the cross 
section per spallation is known then the grammage essentially controls the observed ratio. 
The effects of spallation and energy losses suffered by the secondary nuclei during their 
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residence time in the Galaxy and traversing this grammage, on the average, are to be taken 
into account [13 HH] . 

The models currently in use are exemplified by those of Moskalenko and Strong and of 
Davis et al. [H Ej- The diffusion constant and other transport properties in the models 
are specified in terms of the rigidity (momentum per unit charge) and velocity. These are 
converted to grammage and are shown as a function of kinetic energy, E, per nucleon or 
just kinetic energy for positrons in Fig. 2. These models predict essentially the same energy 
dependence at E > 1 GeV, but differ significantly at lower energies. However, at ~3 GeV 
the grammage of positrons is ~14 g-cm" 2 and ~7 g-cm -2 for Carbon because the Carbon 
nuclei have higher rigidity by a factor of A/Z and a lower velocity compared to the positrons. 
On the other hand, in the alternate model [HI E] discussed here, all the particles are assumed 
to have the same grammage A ~ 1.7 g-cm -2 in the Galaxy. This is also shown in Fig. 2. 
The rest of the energy-dependent grammage left over at lower energies, needed to explain 
the B/C ratio, is attributed to traversal of material in the sources [HJ |9j. 

The following points comparing and contrasting the two models are noteworthy: 

1. The energy dependence of E~ ' 6 of current models [I] will steepen the production 
spectrum S e + ~ e -2 ' 65 to yield a steady state spectrum E~ 3,65 and thence lead to a rapidly 
decreasing positron fraction R e+ ~ E~ x , because the spectrum of the total electronic com- 
ponent is ~ E~ 2,2 at E < 6 GeV and E~ 3A at higher energies. 

2. If we normalize the current models with A e+ (1 GeV < E < 3 GeV) ~ 14 g-cm -2 as 
expected by modeling the B/C ratio, then we over-produce positrons by a factor of ~7 at 
these energies. 

3. This situation is to be contrasted with the model discussed here. The close similarity of 
the spectrum of positrons E e+ (E) and those of the parent primary cosmic-ray nuclei allows 
a good fit to the observations of R e +(E) up to 300 GeV and predicts the expected behavior 
at high energies, reaching an asymptotic value of ~0.7. 

4. Over-production of positrons at low energies is avoided by the fact that the primary 
nuclei have 20 — 30 times higher energy than the positrons they produce, i.e. they will be 
in the range E > 20 — 30 GeV where the residence time in the sources is so short that the 
parent nuclei leak out without significant positron production. 

5. Finally, the ~2 Myr residence time for all cosmic rays at least up to several hundred 
TeV allows one to predict anisotropies consistent with the observations (see [8] for a com- 
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pilation) 5 ~ 3fi;Vn/cn fa 3r/4cr fa 5xl0~ 4 , for the length scale r fa 500 pc. As Strong, 
Moskalenko and Ptuskin have noted ( Fig. 12 in [IS]), the increase of the diffusion constant 
k as E os , which yields A ~ E~ os , concomitantly generates anisotropies that increase with 
increasing energy. This has led to considerable tension between upper bounds on anisotropy 
and that predicted by current models of cosmic ray propagation. 

SUMMARY. 

In summary, we may state that, the production of positrons by nuclear primary cosmic 
rays interacting with the interstellar medium provides a good explanation of the observed 
spectrum and ratio with the total electronic component, provided these particles have a 
residence time of ~2 Myr in the interstellar medium, independent of their energy. The 
prediction of the positron fraction rests on a calculation of energy loss suffered by the 
positrons during their residence for ~2 Myr in the interstellar medium, which leads to 
~ E~ 3S5 for E > 300 GeV. The total electronic component has a spectrum ~ E~ 31 up 
to ~ 1000 GeV and rapidly decreases in intensity at higher energies, so that the positron 
fraction tends to reach an asymptotic value dictated by their production characteristics as 
secondaries, ~0.7 for E ^> 1000 GeV. Because the sources of primary cosmic-ray electrons 
are discrete, their spectrum shows a cut-off at an energy ~1000 GeV, dictated by the distance 
to the nearest source. The spectrum of primary electrons may be understood as the sum of 
the contributions from various discrete sources, with the nearest source dominating at the 
highest energies El EH [20] . 

ACKNOWLEDGEMENT. 

We would like to thank professors M . H. Israel and S. Nussinov who played a key role in 
shaping these comments through discussions of all the key points made in the paper, and 
honing each scientific point to achieve clarity in the presentation. 



[1] M. Aguilar et al, (AMS collaboration), Phys. Rev. Lett. 110, 141102 (2013). 



8 



[2] O. Adriani et al, (PAMELA collaboration), Nature 458, 607 (2009); Phys. Rev. Lett. 105, 

121101 (2010); E. Mocchiitti, 32nd ICRC, Beijing, (2011). 
[3] M. Ackerman et al, (Fermi-LAT collaboration), Phys. Rev. Lett. 108, 011103 (2012). 
[4] L. Moskalenko and A. Strong, Astropbys. J. 493, 693 (1998). 

[5] H. Yuksel, M. D. Kirstler, and T. Stanev, Phys. Rev. Lett. 103, 051101 (2009); S. Profumo, 
Cent. Euro. J. Phys. 10, 1 (2012). 

[6] L. Bergstrom, J. Edsjo, and G. Zaharijas, Phys. Rev. Lett. 103, 031103 (2009). 

[7] M. Ackerman et al, (Fermi-LAT collaboration), Astrophys. J. 761, 91 (2012). 

[8] R. Cowsik and B. Burch, arXiv: 0908.3494vl [astro-ph.COO] 24 Aug 2009. (includes a com- 
pilation of the spectra of electronic component of cosmic rays). 

[9] R. Cowsik and B. Burch, Phys. Rev. D., 82, 023009 (2010). 
[10] D. Gaggero et al, rXiv:1304.6718v2 [astro-ph.HE] (2013) 



[11] P. Mertsch and S. Sarkar, Phys. Rev. Lett. 103, 081104 (2009); arXiv:0905.3152 [astro-ph.HE] 
(2009). 

[12] N. Shaviv, E. Nakar, and T. Piran, Phys. Rev. Lett. 103, 111302 (2009); a.rXw:0905.6904\ 
[astro-ph.HE] (2009). 

[13] B. Katz, K. Blum, J. Morag, and E. Waxman, Mon. Not. Roy. Astron. Soc. 405, 1458 (2010); 
arXiv: 0907.1686 (2010). 



[14] G. Bernard et al, (Fermi-LAT collaboration), arXiv:1207.4670vl [astro-ph.HE] (2012). 

[15] R. Cowsik et al., Phys. Rev. Lett. 17, 1298 (1966). 

[16] R. Cowsik and M. A. Lee, Astrophys. J. 228, 297 (1979). 

[17] A. J. Davis et al, in AIP Conf. Proc. 528, Acceleration and Transport of Energetic Particles 

Observed in the Heliosphere, ed. R. A. Mewaldt et al. (New York: AIP), 421. (2000) 
[18] R. Cowsik et al, Phys. Rev., 158, 1238, (1967). 

[19] A. Strong, I V. Moskalenko, and V. S. Ptuskin, Annu. Rev. Nucl. Part. Sci. 57, 285 (2007). 
[20] J. Nishimura et al, Adv. Space Res., 19, 767 (1997). 



9 



